Fast transfer and efficient coherent separation of bound cluster in the extended 

Hubbard model 
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We study the formation and dynamics of the bound pair (BP) and bound triple (BT) in strongly 
correlated extended Hubbard model for both Bose and Fermi systems. We find that the bandwidths 
of the BP and BT gain significantly when the on-site and nearest-neighbor interaction strengths reach 
the corresponding resonant points. This allows the fast transfer and efficient coherent separation 
of the BP and BT. The exact result shows that the success probability of the coherent separation 
is unity in the optimal system. In Fermi system, this finding can be applied to create distant 
entanglement without the need of temporal control and measurement process. 
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I. INTRODUCTION 

Ultra-cold atoms have turned out to be an ideal play- 
ground for testing few-particle fundamental physics and 
for their increasing technological applications. The ex- 
perimental observation of atomic bound pair (BP) in op- 
tical lattice stimulates many experimental and the- 
oretical investigations in strongly correlated boson sys- 
tems 0-Q~l] • The essential physics of the proposed BP is 
that, the periodic potential suppresses the single particle 
tunneling across the barrier, a process that would lead 
to a decay of the pair. Such kind of BP therefore cannot 
get migration speed through the lattice. In addition, the 
dynamics of the pair separation and formation is also of 
interest in both fundamental and application aspects. 

In this article we study the influence of the nearest- 
neighbor (NN) interaction on the few-boson bound state. 
We will show that, the resonant NN interaction can also 
induce the BP and bound triple (BT) states. These 
bound states are distinct from the known ones, as their 
bandwidths are comparable to that of a single particle. 
This feature allows the full overlap between the scatter- 
ing band of a single particle and the band of a BP or a 
BT, leading to the coherent separation and combination 
processes. We propose relatively simple schemes to im- 
plement such processes. The exact result shows that the 
success probability of the coherent separation is unity for 
a BP while approaches to unity for a BT in the optimal 
systems. Applying the scheme to the extended Fermi- 
Hubbard system, the coherent separation process can be 
utilized to create long-distance entanglement without the 
need of temporal control and measurement process. 



II. FAST BOUND PAIR DYNAMICS 

Let us start by analyzing in detail the two-particle 
problem in an extended Bose-Hubbard model, which is 



simpler than the Fermi one but shares the similar features 
for the issue concerned. The Hamiltonian is written as 
follows: 



H n = -t 
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where a\ is the creation operator of the boson at the ith 
site, the tunneling strength, on-site and NN interactions 
between bosons are denoted by t, U and V. For the sake 
of clarity and simplicity, we only consider odd-site system 
with N = 2Nq + 1, and periodic boundary conditions 
ajv +1 = a\. In this article, we only concern the one- 
dimensional system. Nevertheless the conclusion can be 
extended to high-dimensional system. 

First of all, a state in the two-particle Hilbert space, 
as shown in Ref. [po| . can be written as \ipk) — 

E k ,rf k (r)\^r), With 



i) = -/^E^'^Vac), ( 2 ) 
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here |vac) is the vacuum state for the boson operator 
<Zj. k = 2nn/N, n g [1, AT] denotes the momentum, and 
r e [1, No] is the distance between the two particles. Due 
to the translational symmetry of the present system, the 
Schrodinger equations for f k (r), r g [0, No] are easily 
shown to be 
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where Tf = -2y/2t cos(fc/2) for r = 0, and -2£cos(fc/2) 
for r ^ 0, respectively. Besides, we also have the bound- 
ary conditions f k {-l) = f k (N + 1) = 0. Note that for 
an arbitrary k, the solution of (j4|) is equivalent to that 
of the single-particle Nq + 1-site tight-binding chain sys- 
tem with NN hopping amplitude T k , on-site potentials 
U, V and — 2t cos(fc/2) at Oth, 1th and iVoth sites respec- 
tively. Obviously, in each /c-invariant subspace, there are 
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three types of bound states arising from the on-site po- 
tentials under the following conditions. In the case with 
\U — V\ ^> t, the particle can be localized at either Oth or 
1th site, corresponding to (i) the on-site BP state, or (ii) 
the NN BP state. Interestingly, in the case of resonance 
U = V, and \U\, \V\ t, the particle can be in the 
bonding state (or anti-bonding state) between Oth and 
1th sites, corresponding to a new BP state, called (hi) 
the resonant BP (RBP) state. In this article, hereafter 
we forcus on this type of bound state and refer RBP as 
BP. All the N bound states of (hi), indexed by k, consti- 
tute a bound-pair band. 

In previous works [HI [H| , the bound states of (i) and 
(ii) were well investigated and the corresponding bound- 
pair bandwidths are of t 2 /U and t 2 /V order. Then they 
can be regarded as stationary comparing to the single 
particle in the strongly correlated limit. In order to ana- 
lyze the dynamics of the BP of type (hi), we pursue the 
solution of Eq. (j4j in the case with U = V, and \U\, 
\V\ t, via the Bathe-ansatz method. The BP states 
have the form 



f k (r) * Vk 



-V2e k /(e k -U),(r = 0), 
-sgn^Afc^-^r^ 
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with the spectrum 



e k ~ U ± V2, 
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and 



e k = -2tcos(fc/2) , 



(7) 
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is partial dispersion relation of a single particle. We note 
that neglecting the terms with t 2 /U, the two branches 
of the spectrum combine into a BP band and the corre- 
sponding eigenfunction and energy can be rewritten as 



> = E 
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Observing the above expressions, we hnd that it can 
be regarded as a plane wave of the composite particle, 
BP. Actually, this can be easily understood in terms of 
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FIG. 1: (Color online) Schematic illustration of the coherent 
separation process of a BP in the proposed optimal system. 
It consists of a uniform chain and the embedded resonant 
impurity, (a) An incident BP wavepacket comes from the 
left and scatters off the resonant impurity, (b) The reflected 
waves contain the BP and the single particle wavepackets with 
different speeds. A single particle resides at the impurity. For 
an optimal incident wavepacket, the success probability of the 
coherent separation is unity, (c) Appling the previous scheme 
to a singlet femionic BP, it is found that the separated two 
fermions are maximally entangled. 



the equivalent effective Hamiltonian, which can be ob- 
tained in the quasi-invariant subspace spanned by the 
basis {\l),l € [1,2 AT]} with diagonal energy U. The set 
of basis {\l)} is defined as 



10 = { V* l / 2 J /^l vac ) ' ( even 

a (i-i)/2 a (;+i)/2 l vac ) ' ( odd l ) 



(12) 



The effective Hamiltonian of the ring, restricted to the 
basis {|0} and shifted by U, reads 



2N 

Hbp = ~V2t^2 ( lQ(Z + ll + H.c. 
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which describes a free Bloch particle in a single band. 
However, this BP state is distinct from the previous on- 
site and NN BP because of its large bandwidth, 4^/2t. 
It shows that the speed of a BP wavepacket must be in 
the range of ~ 2-^/2^, which is the reason why we call 
it fast bound pair. This is very significant for quantum 
technologies in the following two aspects: (i) The BP 
wavepacket can be a new candidate as quantum infor- 
mation carrier due to its fast speed, (ii) The bandwidth 
matching between BP and a single particle may lead to 
the coherent separation of the BP. 



III. COHERENT SEPARATION OF BOUND 
PAIR 

Considering a system of ([T]) with the additional chem- 
ical potential /i = U , the scattering band of a single par- 
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ticle fully overlaps the band of a BP ([6]). This fact indi- 
cates that a BP may break in the scattering process by 
a resonant impurity. In order to demonstrate a perfect 
process of the BP separation, we propose an optimal sys- 
tem which is a uniform chain with a resonant impurity 
embedded in. The uniform chain acts as a channel for the 
transport of a BP wave, under the conditions U = V , and 
\U\, \V\ ^> t. The resonant impurity consists of two sites 
with the tunneling strength to , one site of which has far- 
off resonant on-site interaction U s satisfying \U S — U\ S> t 
and another site has resonant chemical potential ji = U. 
This setup is described by the Hamiltonian 



h bps = - ( * 53 °i°<+i + (*° ~ f ) a S ai + Hc - ( 14 ) 

\ i— — OO / 

u„ - u 



n l (—+n i+1 - 1) 



-n (n — 1) + Urix 



To investigate the scattering process of an input BP 
wavepacket from the left, one can establish an equiva- 
lent effective Hamiltonian in the quasi-invariant subspace 
spanned by the following basis {\l),l 6 (—00,00)} defined 
as 
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Eqs. (HZ 
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Acting -ffeps on the subspace, the equivalent effective 
Hamiltonian reads 



"Hbps 



f -1 00 \ 
v (=— 00 i—i/ 



-1)(0| + H.c. 



( 16 ) 

which describes a two connected semi-infinite chains with 
different hopping constants. The on-site interaction U s 
ensures Hbps to be a linear chain, which would be of 
benefit to enhance the success probability of the coherent 
separation. 

The separation process of a BP wavepacket is illus- 
trated in Fig. [1] After scattering, a part of the BP 
wavepacket is reflected by the impurity, while the other 
part of it is separated into two independent particles: 
One of the particle resides at the site with chemical po- 
tential U, while the other particle is reflected to the left 
with higher speed than the reflected BP wavepacket. In 
the aid of the effective Hamiltonian ([T5|. the previous 
process can be reduced to a simple single particle scat- 
tering problem: An incident wavepacket is scattered by 
the joint of the two semi-infinite chains. The reflecting 
wave represents the reflecting BP wave, while the trans- 
mitting wave represents the separated reflecting particle. 
In this sense, the coherent separation probability is equal 
to the transmission coefficient of the effective Hamilto- 
nian, which can be obtained exactly via Green' function 
or Bethe-ansatz method [l3l - [l5| . 

For an incident plane wave of fco, the coherent separa- 
tion probability P (ko) for a chosen system of to = \/2t 



P(ko) 



1 - \/2cos 2 (ko) 



sin fcg v/|cos (2fco) I 



(17) 



for k G (7r/4,7r/2] and P (ko) = for k G (0,tt/4]. 
The sudden death of P (ko) within the region (0,tt/4] is 
due to the mismatch between energies of the two sides of 
the joint. On the other hand, P (ko) also represents the 
particle resident population at the impurity. The profile 
of P (ko) is plotted in Fig. [2] It indicates that P (k ) can 
reach 1.0 at ko = tt/2. 

In practice, this process can be implemented via an 
incident wavepacket instead of a plane wave. Actually, 
the typical wavepacket, a Gaussian wave packet in the 
space can be constructed as 



\Hk ,N c )) 



1 <* 



— ~-(i-N c y+ik i 



10, (18) 



where f2 is the normalization factor. Here N c is the center 
of it in the space {\l}}, ko is the momentum of it. The 
corresponding group velocity in the space {\l)} is v g = 

— 2\/2t sin (ko) which is also plotted in Fig. [5] It shows 
that for the fastest wavepacket with momentum ±7r/2, 
the success probability of the coherent separation can 
approach to unity. In addition, it is shown that such a 
wavepacket is the most robust against spreading [l6[ . We 
define it as the optimal BP wavepacket to demonstrate 
the perfect coherent separation. In the original system, 
it can be written as 



(a]) /V2± la ]a] +1 



|vac) , (19) 



which is obtained from (|18p with wide width (a « 1) 
and ko — ±7r/2 by the mapping rule (fTS")) . Here symbols 
± and h c denote the moving direction (sgn(±7r/2)) and 
the center of the wavepacket. The corresponding single 
particle wave packet has the form 



\<p(±,h c )) 



-(j-h c ) ±if j a i | va( ,^ 



(20) 



Note that |</>(±, h c )) is narrower and slower than that of 
\(f(±, h c )). This is also illustrated in Fig. [T](b) and (c). 

Thus, in the optimal case, the perfect scattering pro- 
cess can be expressed as 



\<f>(+, -00)) |vac) 

r \cf>(-, -00)) |vac) +t\tp(- 



3))4|vac}, (21) 



and the corresponding inverse process as 



4 




0.1 0.2 0.3 0.4 
k (units of it) 

FIG. 2: (Color online) Plots of the separation probability 
P (k) (solid line) and the group velocity v g (dashed line) as 
functions of the momentum of the incident wavepacket. It 
shows that both of them get their maxima at tt/2. 



(a]) 2 /V2|vac) 
a\a) |vac) -> (c| f c^ - c^c^) /y/2 |vac) 
Accordingly we represent the wavepacket as 



(24) 



(25) 



\cj>{±,N c )) -> |0 n (±,iV c )), 
\<p(±,N B )) -► \<p a (±,N c )). 
Nevertheless, from the process 



|0t4.(+, -oo)) |vac) =>■ r \<fo±(-, -oo)) |vac) (26) 

+ 71 O^"' _00 ^ °^ ~ _0 °)) c if) l vac ) 

and the corresponding inverse process 



-oo)) a\ |vac) =*> 
r |</>(-, -oo)) aj |vac) + 1 10(-, -oo)) |vac) , (22) 

where =>■ represents the time evolution. Eqs. (|2ip and 
(f22j) represent the reversibility of the processes, coher- 
ent separation and combination, due to the time-reversal 
symmetry of the system. Here r and t represent the re- 
flection and transmission amplitudes. Then the success 
probability of the coherent separation is |i| 2 . Numeri- 
cal simulation is performed to evaluate the separation 
probability. For an optimal Gaussian wavepacket with 
a = 0.01, we have \t\ = 1.00. Hence we conclude that 
the coherent separation scheme is indeed very efficient 
due to its three advantages, fast transfer, robustness and 
high efficiency. 



IV. CREATION OF DISTANT 
ENTANGLEMENT 

Now we turn to analyze the similar problem in the 
extended Fermi-Hubbard model. The Hamiltonian reads 

/ oo \ 

H F = - \ t ^2 C ia C i+lv + (*0 - t) c\ a C Xa + H.C. 
V i—~oo,a J 

oo 

+ U ^ ( n A n ii + niUi+i) + {U s - U) n Qt n oi + Um, 

%—■ — OO 

(23) 

where c\ a is the creation operator of fermion at site i 
with spin a =t,i and n% — + n^. For similar two- 
particle problem, all the previous conclusions for the Bose 
system are valid for the Fermi system under the following 
mapping rule: 



(|pt(+. -oo)) c f u - \n(+> -oo)) c U) l vac ) =^ 27 ) 

+ r (lVt(-> -°°)) C H ~ _0 °)) c it) l vac ) 

V2t\(t> n (-,-oc)) |vac), 

we can see that the two separated fermions have dis- 
tinct feature from that of bosons, which arises from 
the fact that, state ajaj |vac) is a separable state while 

state (ct^ct | — cj_[ct J |,)/\/2 |vac) is the maximally entan- 
gled state. Utilizing this scheme, one can obtain a pair of 
almost maximally entangled nodes, shared by Alice and 
Bob. This can be used for perfect quantum transport via 
teleportation. 

In addition, when we consider the combination pro- 
cess of two fermions with different spin orientations, it 
becomes a little complicated. We will start our analysis 
from the simplest case in two-particle problem, two paral- 
lel fermions. It is equivalent to the two spinless fermions 
system, in which BP no longer exists. The corresponding 
process are 

|<*t(+, -oo)) c it l vac ) l^t(-) -°°)) 4t l vac ) ' ( 28 ) 
which represents the complete reflection of the incident 
wavepacket. Due to the SU(2) symmetry of the Hamil- 
tonian (|T4"|) . we have 

\<Pt(+, -oo)) I vac) + \ -oo)) c[ t | vac) => 

|</J t (-, -oo)) |vac) + |<^(- -oo)) cj t |vac) . (29) 

Eqs. J26l), J28]) and (Sgj| lead to 

[a \<pt(+, -oo)) + p |^(+, -oo))] c\ f |vac) 

fit 

a |v? t (-, -oo)) c{ t |vac) - -j= -co)) |vac) 

+ 2 (lVt(-> -°°)) c u + I'M - ' ~°°)) c if) l vac ) 

- y (ivt(-) -oo)) c u _ "°°)) c it) l vac ) ■ 

(30) 
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The merit of the above scheme lies in the property of the 
natural time evolution process in an always on system 
without the need of temporal control and measurement 
process. 

V. BOUND TRIPLE 

We now consider the three-particle bound state, BT, 
in the extended Bose-Hubbard model. For the sake of 
simplicity, we directly investigate the system in the limit 
|V| 3> t, \U\, where the perturbation method is applica- 
ble and provides a clear physical picture. The BT state 
is constructed by the three-particle cluster in the con- 
figurations, al_ 1 a\al +1 |vac) and (aj) 2 aj ±1 |vac), which 
possess the same diagonal energy 2V. Since the tran- 
sition strength between them (or the bandwidth of the 
new composite particle) is of the order of t, there may 
exists fast transfer mode in such a system. Accordingly, 
its dynamics obeys the following effective Hamiltonian 



Hbt = -y/2t J2 ( |3i-l)(3i| + |3i)(3i + l| (31) 

i=— oo 

W2|3i + l)<3i + 2| +H.c. 



Resonant Impurity Region 



which is obtained based on the quasi-invariant subspace 
spanned by the basis {\l)} 



a l/3-i a l/3 a l/3+i l vac > ' ( mod ( l > 3 ) = °)' 

10 = \ (4-l)/3j 4+2)/3/^l VaC ) ' ( mod (^ - 1=3) = 0), 

'(/-2)/2 ( a \i+i)/2) /V2|vac), (mod (Z- 2, 3) = 0). 

(32) 

The effective Hamiltonian "Hbt depicts a period— 3 or 
trimerized chain system, which can be further brought 
to a diagonal form 



H BT = ^A fc |fc) {k\, 



(33) 



where k = 2 r Kn/N 1 n G [1, A] is the momentum and 
is the eigenstate of %bt- Here we have used the Fourier 
transformation 



1 N 

X,k) = —= e~ ikl \3l - A), 



(34) 



to construct the secular equation 

-V2t -2te~ lk 
-y/2t -V2t 
-2te lk -V2t 

(35) 

for the eigenstate \k) = 52x=i 7A,fe|A, k). The spectrum 
Afc is solution of A| — 8t 2 Aj, + 8t 3 cos k — 0, which pos- 
sesses gaps at ±7r/3 and ±27r/3 due to the trimerization. 



71, fc \ / 7i,fc 

72, k J = A fc 72, k 

73, k I \ 73, k 
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FIG. 3: (Color online) Schematic illustration of the coher- 
ent separation process of a BT in the proposed optimal sys- 
tem. It consists of two semi-infinite uniform chains with dif- 
ferent chemical potentials, resonant impurity, (a) An incident 
BT wavepacket comes from the left and scatters off the bar- 
rier, (b) The reflected waves contain the BT and the single 
particle wavepackets with different speeds. A single particle 
wavepacket moves along the resonant impurity region, (c) 
Schematic illustration of the equivalent effective Hamiltonian 
governs the evolution. For an optimal incident wavepacket, 
the success probability of the coherent separation approaches 
to unity. 



Nevertheless, in the case of the weak trimerization, the 
eigenstates around ±7r/2 are hardly influenced by the 
trimerization. Then the Gaussian wave packet of the 
form (fT5|) with ko — ±ir/2 should be still fast and ro- 
bust against the spreading. The corresponding optimal 
three-particle wavepacket is denoted as |x(±, N c )). 

The scheme to perform the coherent separation of a 
BT is based on the system depicted by the Hamiltonian 

(0 oo \ 

+2 J2) a ^+ i - ^ ta W + h.c. 
i— — oo i—2/ 

oo oo 

+v rum+i + vY^n,, (36) 

i— — oo i—2 

To investigate the scattering process of an input BT wave 
from the left, one can establish an effective Hamiltonian 
in the quasi- invariant subspace spanned by the following 
basis {\l),l € (—00,00)} defined as 



/ Eqs. m, (I < 2) 

1^- \ oL^Oj-x/Valvac) , (I > 2) ■ [6<) 

The corresponding equivalent effective Hamiltonian can 
be written as 

/ 2 +oo\ 



Hbts = ~V2ti +V2y\ \l-l){l\ -V2t\0){3\ 

\l=-oc 1=4/ 


- (2 - V2) t |3j + l)(3j + 2| + H.c. (38) 



j = -oo 
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which is illustrated in Fig. [3] (c) . The equivalent system 
is two connected semi-infinite chains with a side coupled 
two-site segment at the joint of them. In the optimal 
case, the perfect separation (or combination) process in 
the real space can be expressed as 

|x(+, -oo)) |vac) 

+ t \"") 
r \x(~, -oo)) |vac) + t \<p(+, +00)) al^J |vac) , 

and the corresponding inverse process as 

!</>(-> +°°)) a -i a l |vac) =*> 

r \<p(+, +00)) oLiflto l vac ) + t |x( — , -00)) |vac) . 

This separation process is a little different from that of 
the BP. An incident three-particle wavepacket leaves two 
NN pair at the impurity, while a single-particle remains 
going forward. Numerical simulation is performed to 
evaluate the success probability of the previous scatter- 
ing process. For an optimal Gaussian wavepacket with 
a = 0.01 the separation success probability is 0.97 ap- 
proximately. 



VI. CONCLUSION 



In summary, we studied the BP and BT states in the 
extended Hubbard model. We found that, the resonant 
NN interaction can induce the BP and triple states. They 
are distinct from the known ones in the previous works, 
as their bandwidths are comparable to that of a single 
particle. In other words, the bandwidths of the bound 
clusters can be drastically widened by the NN interac- 
tion. We proposed relatively simple schemes to realize the 
coherent separation of these bound clusters. The exact 
result showed that the success probability of the coherent 
separation is unity for a BP while approached to unity for 
a BT in the optimal systems. We also studied the singlet 
BP in the extended Fermi-Hubbard system. We showed 
that the corresponding coherent separation process can 
be utilized to create long-distance entanglement without 
the need of temporal control and measurement process. 
Finally, we believe that our study will shed more light on 
the future research for multi-particle bound state. 

We acknowledge the support of the CNSF (Grant Nos. 
10874091 and 2006CB921205). 
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